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pointed out that some of their expectations do not survive a more careful examination. An effort 
has been made to make the underlying structure transparent also to those who are not familiar 
with details of loop quantum gravity. 
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I. INTRODUCTION 



The status of loop quantum cosmology (LQC) of FRW models has evolved significantly 
over the last year. Specifically, when there is at least one massless scalar field present, the 
physical sector of the theory was constructed in detail and then used to show that the big 
bang is replaced by a quantum bounce. Furthermore, this singularity resolution does not 
come at the cost of introducing undesirable features such as unphysical matter or large 
quantum effects in physically tame situations. We will begin with a brief summary of these 
results particularly because there is some confusion in the literature on the bounce scenario 
and on the difference between LQC and the Wheeler DeWitt (WDW) theory. 

In these LQC models, the availability of the physical inner product on the space of 
solutions to the Hamiltonian constraint and of a complete set of convenient Dirac observables 
led to a precise construction of suitable semi-classical states. Numerical evolution of these 
states then led to a number of detailed and quantitative results: i) Classical general relativity 
is an excellent approximation to quantum theory until matter density reaches ~ O.Olppi, or, 
scalar curvature reaches ~ — 0.157r/£ P1 ; ii) As curvature increases further, quantum geometry 
becomes dominant, creating an effective repulsive force which rises very quickly, overwhelms 
classical gravitational attraction, and causes a bounce at 1 p ~ 0.41ppi, thereby resolving 
the classical singularity. The repulsive force dies very quickly once the density falls below 
0.41ppi; iii) While the classical evolution breaks down at the singularity, the full quantum 
evolution remains well defined, joining the pre-big-bang branch of the universe to the post- 
big-bang branch through a deterministic evolution; and, iv) Contrary to the earlier belief, 
the so-called 'inverse volume effects' associated with the matter Hamiltonian are inessential 
to the singularity resolution in these models. At the Planck scale, dynamics of semi-classical 
states is dictated, rather, by the quantum geometry effects in the gravitational part of the 
Hamiltonian constraint. (For details, see, e.g., [1-4, 6, 7].) 

These results also addressed two concrete criticisms [8, 9] of the earlier status of LQC. 
First, although it had been demonstrated [10] that the quantum Hamiltonian constraint 
remains well-defined at the putative singularity, in absence of a physical inner product and 
Dirac observables the physical meaning of this singularity resolution had remained unclear. 
One could therefore ask [8]: What is the precise and physically relevant sense in which the 
singularity is resolved? Second, while the early formulations of the quantum Hamiltonian 
constraint [10, 11] were an improvement over the Wheeler-DeWitt (WDW) theory near 
the singularity, they turned out to have a serious flaw: they could lead to a significant 
deviation from classical general relativity even in regimes in which the space-time curvature 
is quite low [2]. In particular, in the k=0 model, more the state is semi-classical, lower was 
the matter density and curvature at which the bounce occurred [2, 3]. In particular, for 
perfectly reasonable semi-classical states, the bounce could occur even at density of water! 
In presence of a cosmo logical constant, serious deviations from classical general relativity 
also occurred in low curvature regions well away from the singularity. These problems can 
also occur without a cosmological constant. As a result, in the k=l model it was then natural 



In early papers (e.g., [1-4]) the bounce was said to occur at p — 0.82ppi. However to calculate the 
quantum Hamiltonian constraint these papers used the lowest non-zero eigenvalue of the area operator for 
the area gap. It was later realized that the corresponding eigenstates are not suitable for homogeneous 
cosmologies. On the space of states which are suitable, the area gap is twice as large, which makes the 
critical density half as large [5] . 
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to ask [9]: does LQC predict that there would be a recollapse from the expanding to the 
contracting phase that can occur at low curvatures in classical general relativity? Absence 
of recollapse would be an unacceptable, qualitative departure from general relativity. 

Because the physical sector of the theory has now been constructed in detail, both these 
issues could be addressed satisfactorily [3, 4, 6]. First, 'singularity resolution' refers to the 
behavior of Dirac observables such as matter density, whence the physical meaning of the 
term is now transparent. Second, not only is there a recollapse in the k=l model but 
there is excellent quantitative agreement with general relativity (for universes which grow to 
macroscopic sizes) [4]. These investigations have also shown that, while the singularity can 
be removed rather easily in LQC, considerable care is needed in the formulation of quantum 
dynamics to ensure that the detailed predictions do not lead to gross departures from classical 
general relativity in tame situations [2, 3]. Careful treatment ensures both a good ultra-violet 
behavior (singularity resolution) and a good infra-red behavior (agreement with general 
relativity at low curvatures) [3, 4]. In particular, in the 'improved' LQC dynamics [3, 4], 
the quantum bounce occurs only when the matter density is ~ 0.41/?pi, irrespective of the 
choice of semi-classical states. Once the curvature is low, evolution is close to that predicted 
by classical general relativity also in the presence of a cosmological constant [3]. 2 

However, even within the confines of these simple models, a number of questions still 
remain. First, numerical simulations and effective equations both imply that the bounce 
occurs when the matter density reaches a critical value p cr i t ~ 0.41p p i. It turned out that 
p cr ; t is insensitive to the sign or value of the cosmological constant and is the same whether we 
consider k=0 or k=±l models [3, 4, 14]. Can one understand the physical origin of p crit l Is 
there an upper bound on matter density in physical quantum states which can be established 
without having to make the assumptions that underlie numerical simulations and effective 
equations? Second, the detailed evolution was restricted to states which are semi-classical 
at late times and was carried out numerically. Therefore, an analytical understanding of 
quantum evolution has been lacking and, in particular, the precise reason for the striking 
differences between LQC and the WDW theory has not been fully understood. While states 
which are semi-classical at late times are the most interesting ones, are they essential for 
the bounce scenario as some authors have suggested? Or, is there a sense in which all states 
undergo a quantum bounce? It has been suggested [15] that the near symmetry (in intrinsic 
time) of volume uncertainties around the bounce point is a consequence of the use of a 
special class of semi-classical states at late times. Is this the case or is it shared by a much 
more general class of states which are (truly) semi-classical at late times? 

Similarly, comparison of LQC and the WDW theory raises a number of questions. There 
is a precise sense [2, 3, 11] in which the quantum Hamiltonian constraint of the WDW theory 
approximates that of LQC. Yet, the physical content of the solutions to the two constraints is 
dramatically different: while none of the states of the WDW theory which are semi-classical 
at late times can escape the big bang singularity, all their analogs in LQC do just that. What 
then is the precise relation between the two sets of dynamics? Can one make statements 



2 These important differences between the improved dynamics (sometimes referred to as the /2 evolution) 
and the older one (referred to as the [i evolution) appear not to have been fully appreciated. Indeed, some 
of the recent discussions [12, 13] portray the two schemes as if they are on equal footing. As discussed in 
section VI, the issue of whether any of the proposed dynamics of LQC can be systematically derived from 
LQG is wide open but it is logically distinct from whether they are viable within the confines of LQC. 
Mixing of these two issues has resulted in some unfortunate confusion in the literature. 
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for more general states? While the geometric part of the WDW Hamiltonian constraint 
is a second order differential operator, that of LQC is a second order difference operator, 
the step size being dictated by the 'area gap' A, i.e., the smallest non-zero eigenvalue of 
the area operator. What happens in the mathematical limit in which A is sent to zero by 
hand? Non-relativistic quantum mechanical systems also admit a 'polymer' representation 
(of the Weyl algebra) which mimics the mathematical structure of loop quantum cosmology 
[16]. For a harmonic oscillator, dynamics in this representation involves a (mathematical) 
discreteness parameter t and is distinct from the standard Schrodinger dynamics. However, 
as was shown in detail in [16-18], it reduces to the standard dynamics in the limit in which 
i goes to zero. Is the situation similar in LQC? If so, what is the a precise sense in which 
the LQC dynamics reduces to that of WDW theory in such a limit? 

The purpose of this (and the accompanying [19]) paper is to address these questions. 
Two key ideas will make this task feasible. First, we introduce a new representation in 
both LQC and the WDW theory in which the operator conjugate to volume is diagonal. 
In this representation, the Hamiltonian constraint of LQC becomes a differential operator 
just as in the WDW theory (for the same reasons as in [16-18]). Second, we will introduce 
a harmonic time coordinate — which is tailored to the scalar field clock — already at the 
classical level. This will simplify the factor ordering of the Hamiltonian constraint relative 
to [3] and make the model exactly soluble. In fact, in the WDW as well as in this soluble 
LQC theory the Hamiltonian constraint reduces just to the 2-dimensional Klein-Gordon 
operator! The physical Hilbert space of both theories is then identical (except for a certain 
global symmetry). How can the theories then lead to strikingly different results on resolution 
of singularity? The answer of course is that the (Dirac) observables are represented by 
distinct operators. Thus, in the new representation, one can compare the two theories 
just by studying the relation between two sets of operators on the same Hilbert space. 
Furthermore, in this representation the expression of Dirac observables is rather simple. 
Therefore, one can analytically compute the expectation values and dispersions. We find 
that in the WDW theory the expectation value of the volume operator on a dense set of 
states goes to zero in the distant past (or future). Thus, for a generic state matter density 
diverges in the distant past (or distant future). In this sense the singularity is unavoidable 
in the WDW theory. In LQC by contrast, on a dense sub-space the expectation value of the 
volume operator has a non-zero minimum and diverges both in the distant past and future. 
Thus, the density remains finite and undergoes a bounce. In this sense the quantum bounce 
is generic and not tied just to semi-classical states. 3 The simplified model also enables us 
to show that matter density has a finite upper bound p sup = \/ / 3/32ir 2 'y 3 G 2 h ~ 0.41p P i 
on the physical Hilbert space (where 7 is the Barbero-Immirzi parameter of LQG). This 
value coincides exactly with that of the critical density p cr i t [3, 4]! Since the bound is 
established analytically and without restriction to states which are semi-classical at late 
times, it provides a deeper understanding of the bounce scenario. 

Finally, we analyze the precise sense in which the WDW theory can be regarded as the 
continuum limit of LQC. We show that given any fixed, finite interval / of internal time we 
can shrink (by hand) the area gap A sufficiently so that the WDW theory agrees with LQC 
to any pre-specified degree of accuracy. However, this approximation is not uniform in /; 
the WDW theory cannot approximate LQC for all times no matter how much we shrink the 
area gap. Furthermore, LQC does not admit a limit at all if the gap is shrunk to zero size. 



3 The accompanying paper [19] analyzes dispersions of the volume operator in this solvable LQC. 
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In this sense it is intrinsically discrete, very different from the 'polymer quantum mechanics' 
of a harmonic oscillator [16-18]. We wish to emphasize, however, that all these results refer 
only to FRW models with a massless scalar field. 

We conclude this section with a general observation on singularity resolution. Recall that 
in classical general relativity, it is non-trivial to obtain necessary and sufficient conditions 
to characterize the occurrence of a singularity (for example, a space-time can be singular 
even if all its curvature invariants vanish). It is therefore not surprising that a generally 
applicable and satisfactory notion of 'singularity resolution' is not available in quantum 
gravity. However, in simple situations such as homogeneous cosmologies the notion of a 
singularity is unambiguous in the classical theory. In the quantum theory of these models 
one can also provide a satisfactory notion of singularity resolution. First, one should have 
available the physical Hilbert space (not just the solutions to constraints) and a complete 
family of Dirac observables, at least some of which diverge at the singularity. Then one 
would say that the singularity is resolved if the expectation values of these observables 
remain finite in the regime in which they become classically singular. In the model studied 
in this paper, we have a well defined physical inner product as well as a complete set of 
Dirac observables both in the WDW theory and LQC In the WDW theory the singularity 
is not resolved because the expectation value of the matter density diverges at the putative 
classical singularities while in LQC the singularity is resolved because the expectation values 
of a complete set of Dirac observables including the matter density remain finite. 4 

The paper is organized as follows. In section II we introduce the new representation in the 
WDW theory. Section III introduces the simplification in the LQC Hamiltonian constraint 
and the new representation in which this constraint becomes a differential operator. Section 
IV compares and contrasts the WDW theory and LQC and section V provides an analytical 
understanding of two striking features of LQC. In section VI we summarize the results and 
compares them with analogous discussions in the literature. 

II. THE WDW THEORY 

It is well known that the WDW theory of the k=0 FRW models with a massless scalar 
field is simple because with an appropriate choice of variables the Hamiltonian constraint 
can be cast as a 2-dimensional Klein-Gordon equation (see, e.g., [2, 20]). The purpose 
of this section is to present this theory in a form that makes its relation to LQC easier 
to analyze. For this, we will need to introduce certain variables and notation which may 
appear contrived from the internal viewpoint of the WDW theory, but which will facilitate 
comparison between the two theories in section IV. 



4 Sometimes apparently weaker notions of singularity resolution are discussed. Consider two examples [13]. 
One may be able to show that the wave function vanishes at points of the classically singular regions of 
the configuration space. However, if the physical inner product is non-local in this configuration space 
— as the group averaging procedure often implies — vanishing of the wave function would not imply that 
the probability of finding the universe at these configurations is zero. The second example is that the 
wave function may become highly non-classical. This by itself would not mean that the singularity is 
avoided unless one can show that the expectation values of a family of Dirac observables which become 
classically singular remain finite there. 
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To construct a Hamiltonian formulation of the k=0 model, one has to introduce a finite 
fiducial cell V and restrict all integrations to it [11, 21]. It is easiest to fix a fiducial flat 
metric q c d and let V be a cube of volume V with respect to q c d- The standard canonically 
conjugate variables are (a, p( a )) for geometry and (4>,p^) for the scalar field. Here a is the 
scale factor relative to q c d- the physical metric q c d being given by q c d = a 2 q c d- To compare 
with LQC, we will make a series of transformations on the geometrical pair (a,p( a )). 

First, it is convenient to fix a fiducial co-triad uj a which is orthonormal with respect to 
q c d and consider co-triads u a which are orthonormal with respect to q c d so that u a = a e u; a , 
where e = 1 if uj l a has the same orientation as the fiducial u a and e = — 1 if the orientation 
is opposite. Second, it is convenient to replace the scale factor a by a variable v which is 
proportional to volume (of the cell V). Let us set 

a 3 V . , 

" = E 2<^ (2 ' 1) 

where 7 is a constant (the Barbero-Immirzi parameter of loop quantum gravity) and = 
(GK) l l 2 is the Planck length. Thus v has the dimension of length and, because of the 
orientation factor e, is not required to be positive; it ranges over (—00,00). A canonically 
conjugate variable is given by 

b = -s^^ (2.2) 
3 V a 2 

which has dimensions of inverse length. In quantum theory, the corresponding operators 
then satisfy the commutation relations: 

[b, z>] = 2% . (2.3) 

To construct the Hilbert space of states, we can either use a representation in which v 
is diagonal or b is diagonal. In both cases the Hamiltonian constraint becomes a linear, 
second order differential operator. Thus both representations are simple and neither has a 
particular advantage over the other. In LQC on the other hand, although the two corre- 
sponding representations are again equivalent, in the v representation physical states turn 
out to have support on a discrete set of values v = 4n\, and the geometric part of the Hamil- 
tonian constraint is a difference operator on states ^(^) [3]. Here n is an integer and A is 
the square-root of the smallest non-zero eigenvalue of area. In the b representation, 
on the other hand, states \I/ (b) have support on the continuous interval (0,7r/A) and the 
geometrical part of the Hamiltonian constraint is represented by a second order differential 
operator. 5 Consequently comparison between LQC and the WDW theory is most direct in 
the b representation. We will therefore let the states be wave functions x(b) also i n the 
WDW theory. (As in [2, 3], symbols with underbars will refer to the WDW theory.) 

In this representation, the Hamiltonian constraint takes the form 

^ x (b,0) = 127rG(b9 b ) 2 X (b,0). (2.4) 



5 An analogy with a particle moving on a circle makes this situation transparent. The b representation in 
LQC is analogous to the ^-representation in which the states ^(0) have support on the interval (0, 2ir) 
while the v representation of LQC is analogous to the representation in which pg is diagonal so the the 
wave- functions ^(pe) have support on discrete values pe = nh 
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In the classical theory, is monotonic along all dynamical trajectories and therefore serves 
as an intrinsic clock. As explained in [2], because of the form of the Hamiltonian constraint, 
this interpretation carries over to quantum theory both in the WDW theory as well as LQC. 
Physical states must satisfy (2.4). Furthermore, they must satisfy a symmetry requirement. 
Under change in orientation of the physical co-triad, (u l a ) — > n(o;*) := — ou l a , the spatial 
metric does not change and, since there are no fermions in the model, physics remains 
unchanged. Therefore, this change represents a large gauge transformation. Recall that in 
Yang-Mills theory physical states belong to irreducible representations of the group of large 
gauge transformations (i.e., to the so-called 'theta' sectors). The same reasoning applies to 
the present case. However, now the large gauge transformation II satisfies II 2 = 1, whence 
wave functions are either symmetric or anti-symmetric under this action. 6 Therefore, it 
suffices to restrict oneself to just the positive (or negative) b-half line. For concreteness, let 
us use the positive half line. 

Then, as in the standard WDW theory one can simplify the Hamiltonian constraint by 
replacing b with y which ranges over (—00,00): 

b 



y := - In — or, equivalently b = b e VT ^ y (2.5) 
Vl2irG b 

where b G is a constant of dimensions of inverse length. The constraint now reduces just to 
the Klein-Gordon equation in (y,4>): 

d\ x(y, <P) = d 2 y x{y, 0) =: -0 x(y, <P)- (2-6) 

The physical Hilbert space can be obtained by the group averaging procedure used in loop 
quantum gravity [2, 22-24]. As usual, the procedure tells us that the physical states can be 
taken to be positive frequency solutions to (2.6), i.e., solutions to 

-^a(^) = Vlx(^). (2.7) 

Thus, if the initial data at the intrinsic time = O is x(y> 0o) = dke~ tky x(k), 
then the physical state is the solution 

X (y,<f>) = -7= / dke- ik y +i ^-^x(k) 
v 2tt 



1 



00 




dk e -^+v) e iHo ^ + / dk e ik(t-y)y e -ik4> 



V2vr J-oo V27T 
=: x L {y + ) + x R (y-), (2.8) 

where y± = <p ± y and the subscripts L and R denote left and right moving states. The 
group averaging procedure also implies that the scalar product is the standard one from 
Klein-Gordon theory: 



(x v x 2 ) P hy = -i / dy x^yA) dj, x 2 (yA) - (fy x^yA)) xJyA) 

J 4>=4> 

POO 

= 2 / dk\k\^(k)x 2 (k). (2.9) 



00 



6 Detailed considerations imply that physical states x(b) used here are Fourier transforms of 3i{v)/v where 
^_{v) are the states used in [2-4]. Because of the extra 1/v factor. x(b) are anti-symmetric in the 
WDW theory (as well as in LQC). But this fact will not play an essential role in the subsequent discussion. 
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We will denote the resulting physical Hilbert space by H^. It is obvious from (2.8) and 
(2.9) that the left and right moving sectors of "Hp^™ are mutually orthogonal. We will see 
in section IV that physically the left moving component \ corresponds to the expanding 
branch of FRW space-times while and the right moving component y , to the contracting 
branch. (Had we worked with the b < half line the correspondence would have reversed.) 
It is useful to note that the inner product can also be written as: 



(X v X 2 ) P hy = 2 dy X^vAo) \id y \ X 2 {vAo) (2.10) 

J — oo 

where the absolute value denotes the positive part of the self-adjoint operator id y . 

Our next task is to define a convenient set of Dirac observables. Since in the classical 
theory V\^ a , the volume (of the fiducial cell) at any instant O of internal time, is a Dirac 
observable, in the quantum theory we wish to define a self-adjoint operator V\^ a on "H™^™. 
With this goal in mind, let us first introduce the Schrodinger Hilbert space "H S ch by freezing 
physical states at = O and evaluating the scalar product (2.10) at = O , and define a 
self-adjoint operator v on it. Since v and b are canonically conjugate, we are led to define 

v := Self-adjoint part of — 2i<9b 

-2z 



= Self-adjoint part of , e~ y/mSv d v . (2.11) 

Using the fact that the operator id y is a positive definite self-adjoint operator on the right 
sector "H^ h of "H sc h an d negative-definite on the left sector "H^. h , one can easily show that for 
all smooth functions x i an d X 2 which fall off sufficiently fast at y — ±oo, we have 

a roc r 

(X V -2i d h x 2 U = - / dy (id yX f) e~^ y (id yX «) ~ {id yX [) e~^ v ^dyX L 2 ) 



+ (idyX*) e-^ y (id yX ^) - {id y x[) e" Vi2 ^ {id y x£ 



(2.12) 



where the integral is performed at = O . By calculating the matrix element of the adjoint 
of — 2i d b and adding we obtain the matrix elements of v: 

(X V *X 2 U = - / — - / dy (id yX ?) e-^y (id yX %) - (id y xt) {id y x\) 
~ v ~ l Vl27rGb J-oo L 

(2.13) 

Therefore, the desired operator v on "H sc h is given by 

v = /L= [P R (e^^y id y )P R + P L {e^y id y )P L ] (2.14) 

b Q v YItxG 

where Pr and Pl are projectors on "H^ h and "H^ h . We will now use this expression to define 
a 1-parameter family of (relational) Dirac observables V\^ a on "Hp^™: 

Thus the action of V\,p is obtained by first freezing the positive frequency solution xiv^) 
at = O , acting on it by the volume operator 1^i\ x \ v\ and evolving the resulting function 
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of y using (2.7). (The operator \u\ is the positive part of v on H sc h-) V\^ a is a well-defined, 
self-adjoint operator on "Hp^T because v enjoys these properties on "H sc h- 

The second Dirac observable is much simpler: the momentum pu) = —ihd^. Since it is 
a constant of motion, it is obvious that it preserves the space of solutions to (2.7) and is 
self-adjoint on "H™^™. Finally, we note that and p^) preserves each of the left and right 
moving sectors. Since these constitute a complete set of Dirac observables on "Hp^™, there 
is superselection and one can analyze physics of each of these sectors separately. 

Let us summarize by focusing on the left moving sector for concreteness. In the more 
intuitive Schrodinger representation, physical states are functions Xl(?/) whose Fourier trans- 
form x L (k) has support on the negative half of the fc-axis and which have finite norm (2.10). 
The matrix elements of the basic operators are given by 



(X L , PWXj)pby = 2h / dy dyX L (yAo)d y x^ L (y,(j) ) and 
(X L , "Uo X^phy = j== - / dyd y x L (yAo) e-^y d y x; L (yAo). (2.16) 

V 127rG D J~oo 

where <f) is an arbitrary instant of internal time. Note that v leaves the left sector invari- 
ant and is positive definite on it. Hence, (x L , 0\cj> o x' L ) P h y = (x L , I^UJ x' L ) P hy Finally, 
dynamics is governed by the Schrodinger equation (2.7). 

Remark: In the b representation, the entire theory can be constructed without any ref- 
erence to the constant b D . In the y representation by contrast the constant b G appears even 
in the final theory through the expression of the operator v. (As we just saw, the inner 
product and the expression of the other Dirac observable p^ is independent of b D also in 
the y- representation.) However, it is easy to verify that change in b G just yields a unitarily 
equivalent theory (as it must). To make this explicit, denote by y' the left side of (2.5) 
obtained by replacing b D with b' . Let us again restrict ourselves to the left moving sector. 
Then the theory based on b Q is mapped to that based on b' Q = ab a by the operator U a : 
U a x(y) — x'(y') '■= x(y' + kxa/y/lSnG). This is an unitary map which preserves p^ and 
sends v to z>', i.e., has the action 



U 



2i 



exp (—Vl2nG y) d v 
VT2^Gb v 



2 1 



U^ 1 = exp (-Vm^Gy') d y , . (2.17) 

V 127TG b 



III. SOLVABLE LQC 

This section is divided into three parts. In the first, for convenience of the reader we 
present a brief summary of LQC including a short explanation of the necessity of 'improved' 
dynamics of [3, 4]. These references began with the Hamiltonian constraint corresponding 
to proper time in the classical theory and changed to the internal time defined by the scalar 
field only after quantization. In the second sub-section we follow an alternate strategy: we 
use a 'harmonic time coordinate' (tailored to the use of the scalar field as an internal clock) 
already in the classical theory. Then the quantum Hamiltonian constraint acquires a slightly 
different factor ordering. But this difference turns out to suffice to enable one to solve the 
quantum theory exactly. (The precise relation between the two strategies is discussed in 
the Appendix A.) In the third subsection we introduce the b representation in which the 
LQC Hamiltonian constraint becomes a differential operator, facilitating comparison with 
the WDW theory in section IV. 
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A. Dynamics of LQC 



As in the WDW theory, because of spatial homogeneity and non-compactness of the spa- 
tial manifold, to construct the Hamiltonian formulation one has to introduce an elementary 
cell V and restrict all integrations to it. (This is necessary also in the path integral treatment 
based on Lagrangians.) Let us then introduce fiducial co-triads and triads which define 
a flat metric q a b and choose V to be a cubical cell whose sides are aligned with the three e". 
Symmetries imply that, by a suitable gauge fixing, the basic canonical pair (A^, E^) of loop 
quantum gravity (LQG) can be chosen to have the form [11] 

A\ = c V~ 1/3 ui, and Ef = p ^/d^ V~ 2/3 e? (3.1) 

where, as before, V is the volume of V with respect to the fiducial metric q a b- The dynamical 
variables are thus just c and p. The factors of V are chosen so that the Poisson brackets 
between them are independent of V: 

{c,p} = ^—- (3.2) 

Note that p ranges over the whole real line, being positive when the physical triad e" has 
the same orientation as the fiducial e° and negative when the orientations are opposite, c 
is dimensionless while p has dimensions of area. The two are related to the scale factor 
via p = e a 2 (where as before e — ±1 is the orientation factor), and c = 7a (on classical 
solutions). 

In the LQC Hilbert space, the eigenbasis \p) of p is adapted to quantum geometry. The 
physical volume of the elementary cell V is simply |p| 3//2 , i.e. V\p) = \p\ 3/l2 \p)- The geo- 
metrical part of the Hamiltonian constraint in the connection variables involves the (density 
weighted) physical triads and the field strength F % ah of the gravitational connection A\. 
The factors involving can be quantized [11] using techniques introduced by Thiemann 
[25, 26] in full LQG. The field strength F x ah is more subtle because a fundamental feature of 
LQG is that while holonomy operators are well-defined, there is no operator directly corre- 
sponding to the connection [27, 28]. Therefore, as in gauge theories, components of F % ah have 
to be recovered by considering holonomies around suitable loops, dividing them by the area 
they enclose, and then shrinking these loops. 

In the earlier LQC treatments [11], this area was calculated using the fiducial metric. 
The geometrical part of the Hamiltonian constraint was then a difference operator with 
uniform steps in p, the step size being dictated by the minimum non-zero area eigenvalue 
2\^3ir'ylp l = A£pj. The resulting dynamics resolves the singularity, replacing the big bang 
by a big bounce. However, detailed investigation [2] showed that this dynamics has several 
unphysical features. First, for semi-classical states, the density at which the bounce occurs 
depends on the expectation value p^) in the state, pbounce = (l/187rG7 2 ) 3 / 2 (l/\/2|p(0)|)- 
Now, is a constant of motion and the higher the value of pu) the more semi-classical the 
state is. (For example, in the closed, k=l model the value of p^ determines the maximum 
volume of the universe V max . The larger the value of p^, the larger is V mauX ). Thus, more 
semi-classical the state, lower is the density at which the bounce occurs, whence the theory 
predicts that there can be gross departures from classical general relativity even at density 
of water! Second, such departures are not confined just to the bounce. In presence of a 
cosmological constant A they occur also in a 'tame' region well away from the classical 
singularity, where a 2 A = \p\ A > 1. Finally, even qualitative features of the final quantum 
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dynamics can depend on the choice of the initial fiducial cell V (see Appendix B.2 in [3]). 
The cell is just an auxiliary structure needed in the mathematical framework. Therefore, 
while it is acceptable that intermediate steps in the calculation depend on this 'gauge choice', 
the final physical predictions must not. 

The 'improved dynamics' of [3, 4] is based on the realization that all these limitations 
stem from the fact that fiducial geometry was used in quantization of F l ah . If area in this 
calculation refers instead to the physical geometry, the subsequent mathematics realigns 
itself just in the right way to cure all these problems. The density p cr i t at the bounce 
is now a constant, p C rit = 

v / 3/(327T 2 7 3 G 2 /l) « 0.41/jpi, which is independent of p w and 
therefore has no implicit dependence on the choice of the fiducial cell V. In presence of a 
cosmological constant, departures from general relativity occur only if the curvature enters 
the Planck regime; general relativity continues to be an excellent approximation away from 
the singularity, improving steadily as the scale factor a increases. Thus, the limitations of 
the older dynamics stem from the fact that it uses a mathematically variable but physically 
incorrect notion of area in defining the operator analog of F l ab . Indeed quanta of area — i.e. 
eigenvalues of the area operator — should refer to the physical geometry not to a kinematical 
background. 

B. The Hamiltonian constraint 

It turns out that the geometrical part of the 'improved' Hamiltonian constraint [3] is a 
difference operator with uniform step size, but now in volume (~ a?) rather than p (~ a 2 ). 
Therefore, it is convenient to use a representation in which states are wave functions ^(^), 
where v is given by (2.1), so that the operator V measuring the volume of the elementary 
cell V is given by 7 

VV(v)=2ir% l y\v\*(v). (3.3) 

As before, the variable b of (2.2) is canonically conjugate to v. However, whereas v is a 
well-defined operator, b is no longer well-defined in LQC because now there is no operator 
corresponding to the connection. Nevertheless, since the holonomies are well-defined, exp «Ab 
is a well-defined operator, with action 

explAb^(r/) = * (i/ + 2A) (3.4) 

where A is an arbitrary parameter with dimensions of length. 

Since we wish to use the scalar field <p as an internal time variable and since <p satisfies 
the wave equation □</> = 0, it is natural to use a harmonic gauge in which the time variable 
r satisfies Or = 0. The corresponding space-time metric takes the form: 

ds 2 = a 6 dr 2 + a 2 df 2 (3.5) 

where a = a(r) is the scale factor. (It is easy to verify that r automatically satisfies Or = 
with respect to this metric.) Because the lapse function is now given by N = a 3 , the 



7 In [2-4], states were chosen to be functions of a dimensionless parameter v in terms of which volume 
eigenvalues are 

(87r 7 /6) 3 / 2 (|v|//04i wh ere K = 2v^/3V3\73. Thus, v is related to the present v by 
v = 2/ v 87T7V3 (v/£pi). This is a trivial change of variables to simplify the expression of the Hamiltonian 
constraint. 
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Hamiltonian constraint becomes 



The key simplification is the absence of inverse powers of the scale factor. 8 This simplification 
is sufficient to make the model analytically solvable. Therefore, we will call the resulting 
model solvable LQC and denote it by sLQC This labeling will serve to distinguish the factor 
ordering used here from that used in [3, 4]. 

The "improved dynamics" procedure used in [3] now leads to the quantum constraint: 

02 t , ,\ o ^ i i sinAb sinAb ~ 

8%V(v,</>) = 3nG\u\—j— H— ^— *M ( 3 - 7 ) 

where A is now the (positive) square-root of the area gap: A 2 = A£ P1 = 2\/37r7^p 1 . Since one 
of the goals of this paper is to spell out the relation between LQC and the WDW theory, 
we want to retain the ability of sending the area gap to zero. Therefore we will leave A as 
a free parameter, A < vA-^pi, rather than fixing it to the LQC value A = vA^pi- 

Next, as in the WDW theory, physical states have to satisfy a symmetry requirement. 
Under orientation reversal II of physical triads e", v — > II(z/) = —v. This is a large gauge 
transformation. As discussed in section II, following the standard procedure from gauge 
theories, physical states belong to irreducible representations of the group of large gauge 
transformations. Since there are no fermions in this model, in LQC \P is assumed to be 
symmetric: II\[ / (z/, <fi) := \&(— v, 0) = ^(v, 0). 

Using this symmetry and writing out the explicit action of operators sin Ab, (3.7) simplifies 

to: 

02 t i i\ o ^ sin Ab sin Ab ~ 



(i/ + 2A)*(i/ + 4A) - 2i/*(z/) + (i/ - 2A)*(i/ - 4A) 



37TG 

= 4P 17 

=: 6 H §(i/,0). (3.8) 

The geometrical part, 0(», of the constraint is a difference operator in steps of 4A. Hence, 
as discussed in [3, 4], there is again superselection: for each e G [0,4A), the space of states 
\I/(z/) with support on points v = e + 4nA is preserved under dynamics. In this paper we will 
focus on the e = 'lattice' which is invariant under II. 



C. Reducing the sLQC constraint to a Klein-Gordon Equation 

For reasons explained in the Introduction, we now wish to work in the b representation 
because the geometrical part of the quantum constraint will also become a differential opera- 
tor. Since ^(u, <p) have support on the 'lattice' v = 4nA, and since b is canonically conjugate 



This strategy was not used in [3, 4] because this choice, N = a 3 does not have a natural analog outside 
homogeneous cosmologies while the choice N = 1 used there (which corresponds to using of proper time) 
is generally available. 
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to v, their Fourier transforms \?(b, </>) have support on the continuous interval (0,7r/A): 



tf(b,0):= V e^ b ^(z/,0); so that = - / 7 db e~^ b tf(b, 0) . (3.9) 



From the form (3.7) of the constraint it is obvious that it would be a second order differential 
operator in the b-representation. To facilitate comparison with the WDW theory, let us set 
x{y, 0) = (A/7rz/)\I/(i/, 0). Then, on x(b,0), the constraint (3.7) becomes 



/ sin Ab 



9jx(t>,$ = 127rG 1—^-3,1 X (b,0) (3.10) 

which is strikingly similar to the WDW equation (2.4). Note however, that we did not arrive 
at (3.10) simply by replacing b in the expression of the classical constraint by sinAb/A as 
is often done (see, e.g., [15, 30]). Rather, (3.10) results directly from the 'improved' LQC 
constraint if one begins with a harmonic time coordinate already in the classical theory. 
To simplify the constraint further, let us set 

x = - 1 m(tan — ), or b = - tan" 1 ( e ^^ x ) (3.11) 

so x ranges from — oo to oo. Then (3.8) becomes just the Klein-Gordon equation 

d\ X (x, ct>) = d 2 x X (x, <P) =: -6 X (x, <P) ■ (3.12) 

Therefore we can repeat the discussion of section II. The physical Hilbert space is again 
given by positive frequency solutions to (3.12), i.e. satisfy 



-%(^i) = v6x(j;i). (3.13) 

We can again express the solutions in terms of their initial data and decompose them into 
left and right moving modes x( x i 0) = Xl( x +) + Xr( x ~) ( see Eq (2.8)). The physical inner 
product is given by (2.9) or, equivalently, (2.10). Thus, the physical Hilbert space of LQC 
is exactly the same as in the WDW theory and indeed the action of the operator on 

physical states is the same: p^X — —i^d^x = \f~^tX- 

How can there be difference in physical predictions, then? Recall that to obtain phys- 
ical predictions, we need the action of a complete set of physical observables and, as we 
will now show, the volume observable V\$ in LQC has a different form from that in the 
WDW theory. Consequently, although the Hilbert spaces are mathematically the same, the 
physical meaning of a given positive frequency solution is different in the two theories. 

Let us begin with the operator v on the Schrodinger Hilbert space % sc h spanned by the 
initial data x( x i^o) t° (3.13), where the inner product is obtained by evaluating (2.10) 
at <f) = (j> . As in the WDW theory, v is again the self-adjoint part of —2x8^. However, 
while —2idi, = (— 2z/ v / 127rc/b ) exp(— Vl27r Gy) d y in the WDW theory, it is now given by 
{—2i\ / yY2nG) cosh(A/l27rGx) d x . Therefore, repeating the procedure we followed in sec- 
tion II, on % sc h the operator v is given by: 

2A 

v = [P /? (cosh(v / 127rGx) id x ) P R + P L (cosh(Vl2nG x) id x ) P L ] . (3.14) 

V \2ixG 



13 



Therefore, as in the WDW theory, the corresponding volume operator on "H p h y is given by 

VU xM)=e* /S <-+-*° ) (27r T 4 1 |i>|) xMo) ■ (3.15) 

A second difference from the WDW theory arises because of symmetry conditions on 
the physical states. Recall that to qualify as a physical state ^(u, 0) has not only to be a 
positive frequency solution but must also satisfy v, 0) = ^(u, 0). This translates to the 
condition x(— x, 0) = — x(a;,0). Therefore x{ x i 4>) has the form 

x(x,0) = -^=(F(x + )-F(x_)) (3.16) 

for some function F (such that F{x±) are positive frequency solutions to (3.12)). Conse- 
quently, in contrast to the WDW theory, the right and left sectors are not superselected. 
However, since the full information in any physical state xi x i 0) is contained in F, we can 
conveniently describe the "H p h y in terms of positive frequency, left moving solutions F(x + ) 
(or, right moving solutions F(x_)) alone, which are free from any symmetry requirement. 
In this description, the scalar product is given simply by 

/oo 
dx F 1 (x + )d x F 2 (x+) (3.17) 
-oo 

where the integral is evaluated at a constant value of 0. 

Let us summarize the basic structure using a Schrodinger representation which will be 
useful in the next section. For any given A, the LQC Hilbert space can be taken to be the 
space of functions F(x) whose Fourier transform F(k) has support only on the positive half 
line and whose norm given by (3.17) is finite. Matrix elements of the basic observables pu,\ 
and v are given by: 

/oo 
dxd x F 1 (x + ) d x F 2 (x + ) and 
■oo 

4A f°° 

(Fi, % F 2 ) phy = J dxd x F 1 {x + ) cosh(Vl27rG'x) d x F 2 {x+) (3.18) 

where the integral is evaluated at = O . The states evolve via the Schrodinger equation 

-id 4> F(x + ) = VOF(x + ) (3.19) 

so that, as in the WDW theory the effective Hamiltonian is non-local in x. Thus, the only 
difference from the WDW theory lies in the expression of v and hence of the volume operator. 

Remark: A quick way to arrive at the constraint (3.10) in the b-representation is to 
write the classical constraint in terms of the canonical pair (b, v) and then simply replace 
b by (sin Ab)/A and v by —2id b . While this so-called 'polymerization method' [18, 30] 
yields the correct final result, it is not directly related to procedures used in LQG. Since 
v is a geometrical variable, its quantization could be carried out using quantum geometry 
techniques. However, b has no natural analog in LQG. In particular, since it is not a 
connection component, a priori it is not clear why the wave functions have to be even 
almost periodic in b nor why A should be related to the area gap. For a plausible relation to 
LQG one has to start with the canonical pair (A z a , E" 1 ), i.e., (c,p), mimic the procedure used 
in LQG as much as possible, e.g., along the lines of [3], and then pass to the b representation 
as was done here. Once this is done, a posteriori it is possible, and indeed often very useful, 
to use shortcuts such as b — > (sinAb)/A. 
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IV. THE WDW THEORY AND LQC: SIMILARITIES AND DIFFERENCES 



This section is divided into two parts. In the first we show that there is a precise sense in 
which the singularity is generic in the WDW theory while it is generically replaced by a quan- 
tum bounce in sLQC In the second we spell out a precise sense in which the WDW theory 
approximates sLQC and the sense in which this approximation fails. 



A. Singularity versus the quantum bounce 

In spatially homogeneous, isotropic space-times the only curvature invariant is the space- 
time scalar curvature which is proportional to the matter density. Classical singularity is 
characterized by the divergence of these quantities. Now, the matter density is given by 
p = p^/2V 2 and since there is no potential for the scalar field, is a constant of motion 
in our model (also in the quantum theory). Therefore, to answer the question of whether 
a given quantum state is singular or not, one can calculate the expectation values of the 
Dirac observable V<j,. If (V"^) goes to zero (possibly in the limit as <f> tends to ±00 as in 
classical general relativity), one can say that the quantum state leads to a singularity. If 
(V\,p) diverges in both directions, i.e., as <p ~ > ±00, but attains a non-zero minimum, one 
can say that the state exhibits a quantum-bounce. While this is a rather weak criterion, it 
is applicable to any state (in the domain of the volume operator). In this sense, it enables 
us to test the generic behavior in the WDW theory and in sLQC, thereby going beyond the 
analysis of [1-3] which was restricted to states which are semi-classical at late times. 

Let us begin with the WDW theory. As noted in section II the left and right moving 
sector are superselected and can be analyzed separately. For definiteness we will focus on 
the left moving sector. Then, for any x{v+) we have 



VV2irGb 
2 

pi 



dy 



dy 



Vl2nGb 



=: V a e 



Vl2nG(t> 



(4.1) 



where V is a constant determined by the solution and can be calculated from initial data 
at any instant of time. Thus, for any state x L (y+) m the domain of the volume operator, 

the expectation value {V\$) tends to 00 as the internal time tends to 00 and goes to 
zero as 4> tends to —00. In this sense the left moving sector corresponds to the expanding 
universes. It is obvious from the above calculation that on the right moving sector x R (U-) 
the situation would be reversed whence it represents contracting universes. In the first case 
there is a big bang singularity and in the second, a big crunch singularity. Indeed, states 
which are semi-classical in an epoch when the density is very low compared to the Planck 
density are known to follow the classical trajectories into either the big bang or the big- 
crunch singularity [3]. We have shown that this qualitative behavior is generic: for a dense 
subset of states, expectation values (V"|^) of the volume operator evolve into a big bang or 
a big crunch singularity. This calculation also shows that the matter density, defined as 



15 



P '■= (P(</))) 2 /2(V r ) 2 diverges as <fi — > — oo (resp. </> — > oo) in the left (resp. right) moving 
sector. Furthermore, the same considerations apply to the expectation values of the density 
operator /^because, as shown in section VA, its expectation value has the same behavior 
as p. 

Let us carry out the same calculation in sLQC. Consider any state x{ x i 0) = 
(1/V2) (F(x+) - F(x_)). Then, we have: 



(x, vUx) 



phy 



2tt74i(x, I^UoX) 



phy 



Vl2nG 
87T74iA 

which can be written as 



dx [d x xi cosh(v / 127rGx) d x \ L + 8 x xr cosh(v / 127rGx) d x xn] 
cosh ( V YI-kG (x. 



dx 



dF 



dx_i 



(x,v\<i>x) 



phy 



Vl2wG<f> 



+ V.e 



<t>)) 



where 



V + 



POO 

/ dx+ 


dF 


2 -1 

e q=Vl27rGx+ 


_ J —00 


dx + 





(4.2) 



(4.3) 



are again constants associated with the solution which can be calculated from the initial 
data at any fixed (f). Note that V± are strictly positive (compare [15]). In contrast to the 
WDW theory, (V^) tends to infinity both in the distant future and in the distant past. 
Furthermore, it has a unique local minimum which is therefore also a global minimum. In 
this sense the bounce picture of [3, 4] is robust. It is not restricted to states which are 
semi-classical at late times but holds for all states in the domain of the volume operator. 
Similarly, since the result is obtained analytically, an apparent concern (which arose from 
some observations made in [31]) that the presence of the bounce may be related to intricacies 
of numerics is ill founded. The bounce point and the minimum volume are given by: 



1 , V. 
— , log — , 
2 V / 127G S V 



and Vrr 
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(4.4) 



where is the norm of the state x- 
Next, we note that 

(VU)^ = V min cosh(v / 127G0) . (4.5) 

Therefore, the internal-time evolution of volume is exactly symmetric about the bounce 
point for all states. For reasons just discussed one would expect that the symmetry would 
hold also in LQC to an excellent degree of approximation. This symmetry was already ob- 
served in [3] through numerical and effective equation methods and in [15] using a simplified 
model. However, all these analyses were tied to various notions of semi-classicality. Our 
results overcome the concern that the symmetry may not extend to general states. Our 
considerations here refer only to expectation values. However, the same techniques are used 
in [19] to calculate fluctuations of the volume operator and analyze its properties. 

Remark: In the WDW theory the left and right sectors are superselected and should 
therefore be considered separately. Nonetheless, since in sLQC all physical states x satisfy 
\/2x(x, 0) = F(x + )— F(xJ) for some F, one may be tempted to consider states x(y, 0) which 
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are similarly anti-symmetric in y also in the WDW theory. This strategy would however be 
physically incorrect because of two reasons. First, the anti-symmetry condition in x arises in 
sLQC because physics is symmetric under reflections of spatial triads Ef. This condition has 
already been incorporated in the WDW theory in arriving at the y representation. Secondly, 
if one were nonetheless to impose anti-symmetry in the WDW theory, one would find that 
this sector does not have any semi-classical states at very late (or very early) times because, 
in this sector, states which are peaked at large volume would have very large — rather than 
small — extrinsic curvature. 

B. Approximating sLQC with the WDW theory 

There is a precise sense in which the WDW equation approximates the LQC evolution 
equation [2, 3, 11]. However, as is well known, the issue of whether solutions to two such 
equations approximate each other is logically quite distinct. In particular the equations may 
be such that initial data satisfying the assumptions needed to show that the two equations 
approximate each other may evolve out of this regime. If this occurs, solutions of the two 
theories may be close to each other initially but could eventually differ drastically. In this 
section we will show that this is precisely what happens in our case. 

At first this may seem perplexing because in sections II and III we showed that the 
physical states in both theories satisfy the (positive frequency) Klein-Gordon equation. Isn't 
then the evolution identical? However, as pointed out in section III C that in quantum 
theory the physical content of a state lies not in its functional form but in its interplay with 
observables. While the action of p^ is the same in the two theories, the action of is 
distinct (compare (2.16) and (3.18)). Therefore the theories are inequivalent. The question 
for us is: Can one arrive at the WDW theory by shrinking the area gap of sLQC, i.e., by 
letting A go to zero? Again, one cannot naively set A to zero: A appears linearly in the 
expression of the volume operator (see (3.18)) and setting it to zero will make the operator 
itself zero. Rather, we have to work with states and operators together. 

Let us consider the 1-parameter family of Dirac observables V\ ( f > . A necessary condition 
for the WDW theory to approximate sLQC is that, for any state of the WDW theory, there 
should exist a state of sLQC(a) such that the expectation values of volume operators in the 
two states can be made arbitrarily close to one another by choosing a sufficiently small A. 9 
More precisely, given e > 0, there should exist a 5 > such that for any state xAVi^) °f 
the WDW theory, there is a state X w ( x ^ 4>) °f sLQC(a) such that 

\(VU) ( x)-(VU) wdw \ <e (4.6) 
for all A < 5. Now, the forms (4.1) and (4.2) of the two expectation values 

(^U)wdw = Voe^^ and (V^) (A) = V+e^ 5 * + V_ e -^^* (4.7) 
immediately implies that if (4.6) is to hold for all positive then 

V = V+ and V_ < e (4.8) 



9 In the k=0, spatially non-compact context considered here, one must introduce a finite cell to speak of 
volume. However, ratios of the volume expectation values (V - ^) at two different times or, at the same 
time but computed using sLQC and the WDW theory are independent of the choice of the cell. 
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where V is defined by (4.1) using x L {Vi 0) an< ^ V± by (4.3) using x (A) (x,0). However, since 

V_ is necessarily positive, it follows that |(V"|</,)(a) — (V^)wdw| will grow unboundedly as 
0—7- — oo no matter how small 5 is. Thus, the necessary condition for the WDW theory to 
approximate sLQC is violated: the global dynamics of the two theories is very different. 

While the simple argument above establishes the main result of this section, it is 
instructive to examine some details by explicitly constructing a map from XriVifi) to 
X w (x,<f>) = (1/ V2)(F^x){x + ) — F^){ x -))- For this it is convenient to work in the Schrodinger 
representation at time, say, = 0. Then our task is to construct from the initial 

data 0). Fix an e > and consider the family of Schrodinger states X L (y) °f the 

WDW theory for which V_ of (4.3) (obtained by using x L f° r F) satisfies V- < e. Then, 
considerations in the remark at the end of section II suggest that we set 

1 2 

F(\) (x) := X T (x + fi) where fi = -j== In — . (4.9) 

Then for any A, the pairs x(u) an d F^(x) are Schrodinger states which are 'close' to one 
another at = in the sense that: 

i) They have the same norm: \\x{y)\\p$? = \\F(x)(x)\\p hy ; and, 

ii) They have the same expectation values for any power of p^: (pj^) d = {P(<p)), x y an d> 

iii) Their volume expectation values are close: | (V\ ( p) wdw — (V"|</>)(a)I < e. 

Note that since \x depends on A, for each A the prescription chooses a different initial state 
in sLQC(a). If x(y) were peaked at some value y , each Fn\(s) will also be peaked at some 
x Q . However, as A decreases this peak would shift progressively to x = —oo. This change 
of the functional form of Fn) (x) with A is the 'renormalization flow' required to ensure that 
physics of the initial state remains the same at each scale A. 

Let us now evolve these states. Since is a constant of motion, the non-trivial dynamics 
is in the expectation values of the volume operator. (4.7) immediately implies that while 
the difference between the two sets of expectation values will in fact shrink as <fi increases, in 
the distant past it will grow unboundedly! Thus, the two sets of dynamics approximate each 
other well only on a semi infinite time interval (which is finite in the negative direction). 

What happens when we shrink A? Since 



V- 



47T7giA 
v / 127G 



/■oo 

/ dx 


dF 


2 -i 

e Vl2nGx 


. J — oo 


dx 





(4.10) 



and since F^)(x) = x( x + A 4 ) where \i monotonically decreases and /x — > — oo as A — > 
0, V- decreases steadily in this limit. Therefore, it follows from (4.2) that the condition 
|(^U)wdw(0) — (^U)(A)(0)| < e is satisfied further and further in the distant past. Therefore, 
given an e > and any semi-infinite time interval I = (— O , oo), however large <p may be, 
there exists a A > such that the dynamical evolution of (V) in sLQC(a) dynamics remains 
within e of that in the WDW theory on I. In particular, as observed in [3], as we shrink A, 
for the same initial data at = 0, the bounce time B is pushed further and further into 
the past. 

To summarize, the LQC dynamics does not reduce to the WDW dynamics as the area 
gap shrinks to zero. However, if one is interested only in semi-infinite intervals such as 
(— O , oo), one can recover the WDW dynamics on choosing by hand a sufficiently small A. 
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V. TWO STRIKING FEATURES OF SLQC 



This section is divided into two parts. In the first, we show that matter density admits 
an absolute upper bound p sup on the physical Hilbert space which, furthermore, equals the 
critical density p crit found in [3, 4] using effective equations and numerical evolutions of the 
exact equations. In the second, we show that there is a precise sense in which sLQC is a 
fundamentally discrete theory. 



A. An absolute upper bound on matter density 



The fact that every LQC physical state undergoes a bounce at a positive value of (V"]^) 
provides a sense in which the singularity is resolved in LQC. We will now show that there is 
a much stronger sense in which the resolution occurs: the spectrum of the density operator 
p\(j, is bounded above on the full physical Hilbert space "H p hy Note that the boundedness of 
volume by itself is not sufficient to imply the boundedness of p]^ because the operator p^) 
does not have a finite upper bound on "H p h y . 

In the classical theory, the scalar field density is given by p\^ = (l/2)(p^/V\ ( p) 2 . Since 
the operators p^ and V\ ( j > do not commute in quantum theory, it is natural to define the 
density operator as 



P\4> 



1 A\ 2 

2 l * 



where A\^ = (V^y 1 ' 2 p^iV^y 1 ' 2 



(5.1) 



Now, because p^ is a constant of motion, in any given physical state one would expect 
the time-dependent observable p\^ to reach its maximum value p B at internal time <fi = 4>b- 
Numerical simulations have shown [3, 4] that for states which are semi-classical at late times 
the value of p B is remarkably robust: p B ~ 0.41p P i for all the states considered in models 
with and without a cosmological constant and also for the k = ±1 cosmologies. (It was 
therefore called critical density and denoted p C rit)- It is natural to ask if this fact has an 
analytical explanation in sLQC. We will now show that the answer is in the affirmative. 

Let us first compute the expectation values of the operator A using those of p^ and 
given in (3.18) 
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where 



xM) = v 



V2 



(F(x+) - 



(5.2) 



(5.3) 



and the integrals are performed at = < 
of the the two integrals is bounded above by 1 
in LQC, we obtain 



(A>< 



Since cosh(v / 127rGa;) > 1, it follows that the ratio 



Therefore, using the fact that A 2 = 2\^3 i y£ 



pi 



47T7 2 G 



1 

A 



(5.4) 
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This implies that the spectrum of is bounded above by the right side of (5.4). Therefore 
the spectrum of p\^ = A 2 \ ( p is also bounded by 

Psup = 8^G h = 32^ W °- 41PP1 ' (5 - 5) 
where in the last step we have used the value 7 ~ 0.24 for the Barbero-Immirzi parameter 
that led to p cr ;t ~ 0.41p p i in [3, 4, 14]. We wish to emphasize that this is an absolute bound 
on the spectrum of p\$ on the entire physical Hilbert space "H p h y ; there is no restriction that 
the states be, e.g., semi-classical. Note also that a factor of 10 in the value of 7 would change 
p sup (and Pcrit) by a factor of 10 3 . The fact that the value of 7 obtained from the entropy 
calculation yields p sup ~ ppi points to an overall coherence of LQG. 

Finally, we could also have defined, as a measure of mean density, the quantity p = 
(P(<A)) 2 /2(^U) 2 - It is easy to verify using the reasoning that led us to (5.5) that p is also 
bounded above by p sup . 

Two questions arise immediately: Is this supremum attained? And, how does it relate to 
the upper bound p B on density along individual dynamical trajectories! The fact that the 
value of p sup is the same as p crit suggests that the answers to the two questions are related. 
This is indeed the case. 

Let F(k) be a smooth function satisfying the following conditions: 

F(k) = l e - ^ (fc 2 fco)2 e fe ° if k > eVG > and F(k) = if k < , (5.6) 
k 

where e <C 1 <C k /\/~G. This is a semi-classical initial state peaked at k = k a and x = x a . 
Let F(x) be its Fourier transform. Then x{ x i 4>) = 0-/V%) (F(x + ) — F(x-)) is a physical 
state. One can readily calculate the expectation values and {V\^} in this state. Since 

(j>u)) is a constant of motion and (V\^) attains its minimum, V B , at the bounce point, as 
expected p attains its maximum also at the bounce point. It is given by: 



Pb ~ Psup 



1-0 



Gh 2 



pU + ( A P(<«) 2 



(5.7) 



where we have used the fact that Ap^ = h/3. 

Thus, for semi-classical states p B is very close to p sup - Typical numerical simulations of 
[3] used such semi-classical states with (p(^)) = 5 x 10 3 h (in the classical units G — c — 1) 
with the relative dispersion Ap^/ (p^) of 2.5%. In this case, the above calculation shows 
that p sup — p B = O(10~ 4 ), whence the value p cr i t = 0.41ppi of that numerical simulation is 
consistent with the above analytical calculation in the simplified model. Note incidentally 
that even this state represents a universe that is very quantum mechanical. In the k = 1 case 
for example, a state with these parameters represents a universe which grows to a maximum 
radius of only ~ 25£p\ before undergoing a classical recollapse. If we use values of {p(<f>)) and 
Ap(0) that correspond to universes that grow, say a megaparsec size in the k — 1 case [4], 
Psup would agree to p B to 1 part in 10 230 . These considerations show that on % p h y , p B can 
come arbitrarily close to p SU p- 



B. Fundamental discreteness of sLQC 

We saw in section IV B that the WDW dynamics does not result in the limit A — > of 
sLQC(A). A natural question then is whether sLQC(a) admits any limit at all as the area 
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gap shrinks to zero. We will now show that the answer is in the negative. The underlying 
conceptual idea is simple: the one-parameter family of theories sLQC(a) will admit a limit 
as A — > if, given any e > there exists a 5 > such that the predictions of all sLQC(a) 
with A < 5 agree within e for all times <fi. We will argue by contradiction. 

Suppose then that the limit exists. Fix an e and let A, A' < 5. Fix a normalized state % 
in sLQC(a). Then by assumption sLQC(a') admits a normalized state x , which is 'close' 
to it. In particular then, 

|<^>(A)-(^>(A')l<e, (5-8) 
and, to keep transparency of equations we will assume that (p(<^))(A) — (P(4>))(y)- W Now, 

S1I1CG 

(VV ( A) = ^ + (A)e^ /T ^* + V„(A)e" v/T ^^, (5.9) 

Eq (5.8) can be satisfied for all (f> if and only if the V±(A) = V±(X'). It then follows from 
(4.4) that the densities at the bounce points of the two theories are equal: 

Pb(A) " 8 V + (X)V.(X) ~ 8 V + (\')V-(\') ~ Pb(A } • (5 ' 10) 

Now fix A' < 5, choose A = X'/N for some N > 1 and let x (A) be a semi-classical state with 
Pb(X) = Psup(A) — e. As we saw in section VA, we can choose % so that e is arbitrarily 
small. But then (5.10) and (5.5) imply 

Pb(X') = p B (X) = p sup (A) - e = N 2 Psup (X') - i > p sup (A') (5.11) 

if we choose to be sufficiently large. But this is impossible since p S up(A') is an absolute 
upper bound on density in the physical Hilbert space of the sLQC(a') theory. This implies 
that our initial assumption that the limit exists is invalid. sLQC does not admit a continuum 
limit; it is a fundamentally discrete theory. This establishes the main result of this section. 

It is useful to note that, by replacing b G A with A/A' in (4.9), the procedure introduced in 
section IV B leads to a rather natural 'renormalization flow' relating states of the 1-parameter 
family of sLQC(a) theories. Under this flow, if the initial x (x) { x ) is semi-classical in sLQC(a), 
the state X {x/) { x ) w iU a l so be semi-classical in sLQC(y). The classical trajectories on which 
the two states are peaked in the future will be very close to one another. Each state will also 
remain peaked on a classical trajectory in the (pre-bounce) distant past. However, these 
two classical trajectories will be very different and diverge away from each other in the past 
evolution. Thus, while there is a bounce in each sLQC(a), the pre-bounce dynamics is vastly 
different for different values of A. 

The non-convergence of sLQC in the A — > limit shows that sLQC is qualitatively dif- 
ferent from polymer quantum mechanics [16-18] and lattice gauge theories. The polymer 
particle example was introduced as a toy model to probe certain mathematical and con- 
ceptual issues and by itself does not have direct physical significance. While it has certain 
similarities with LQC, there are also important differences. Because there is no positive 



It is completely straightforward to extend the argument to allow these expectation values to be different 
but within e of each other. Requiring an exact equality is also motivated by the fact that p<^\ is a 
constant of motion and, even within a single sLQC(a) theory, semi-classical states with different values of 
depart from each other unboundedly under evolution. 
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frequency requirement in the polymer particle theory, the physical inner product and the 
Hamiltonian are local. Arguments that are successful in that example simply break down 
because of non-locality of sLQC. Lattice gauge theories are meant to be controlled approx- 
imations to the continuum theory. As in numerical analysis of PDEs, the discreteness is 
introduced just as an intermediate mathematical simplification. Therefore, absence of a 
continuum limit would be tantamount to non-viability of that theory. In LQG on the other 
hand the continuum is only an approximation. Not only is there no a priori reason for the 
theory to make sense as the area gap is shrunk to zero, but within the LQG framework it 
would be hard to make physical sense of the limit if it existed. Discreteness at the Planck 
scale is a fundamental and essential ingredient of the theory. The continuum emerges on 
coarse graining; ignoring the fine structure of quantum geometry because one is not inter- 
ested in phenomena at the Planck scale is very different from taking the naive continuum 
limit A —> which corresponds to washing out quantum geometry at all scales. 

VI. DISCUSSION 

Detailed analysis of the physical sector of LQC has revealed that the older quantum 
Hamiltonian constraints [10, 11] have serious drawbacks already in the FRW models coupled 
to at least one massless scalar field [2]. In particular, they lead to physically unacceptable 
breakdown of general relativity in completely 'tame' regimes. The 'improved' constraint 
operator of [3] is free of these drawbacks. Within the FRW models, this improvement is 
robust in the sense that it extends to situations in which there is a non-zero cosmological 
constant and/or non-zero spatial curvature corresponding to the k=l cosmologies [3, 4]. 
Recent analysis has shown that its bounce picture is also robust with respect to the inclusion 
of a phenomenologically viable inflationary potential [32]. Thus, the suggestion (see, e.g., 
[13, 15]) that the bounce would not persist once a potential is included has turned out to 
be incorrect. 

However, so far predictions from 'improved' dynamics were all made numerically using 
states which are semi-classical at late times. Three somewhat different classes of semi- 
classical states were used and numerical simulations were carried out in two distinct ways, 
using a variety of values of the parameters involved. While these features had established a 
certain degree of robustness of results, in absence of an analytical understanding of dynamics, 
it was not obvious whether the results would continue to hold for generic states. Indeed, it 
has been suggested that they may not [15] and concerns have been expressed that numerical 
subtleties may affect the robustness of the bounce [31]. 

In this paper we obtained an analytically soluble model — sLQC — by adapting the theory 
to the scalar field clock already at the classical level. This allows one to consider generic 
states and analyze physics without recourse to any numerics. The question then is: Does 
the quantum bounce persist generically and, if so, does it continue to be approximately 
symmetric about the bounce point? Or, are these features restricted just to semi-classical 
states? A natural avenue to explore these questions is through dynamics of the expectation 
values of the volume operator. We found that these undergo a quantum bounce for all 
states and furthermore the bounce continues to be symmetric. Thus the bounce and its 
salient qualitative feature are quite robust in the cosmological model under consideration. 
More importantly, we could show that matter density has an absolute supremum p sup on 
the physical Hilbert space, given by: p sup = \/3/32ir 2 'y 3 G 2 h ~ 0.41p P1 . This is precisely the 
value of the critical density p crit [3, 4] at which the bounce occurs in numerical simulations 
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and in solutions to effective equations, both of which are, however, based on states that 
are semi-classical at late times. Thus the raison d'etre and the observed robustness of the 
somewhat mysterious p crA was clarified analytically. 

We could also use the model to analyze the relation between the WDW theory and sLQC. 
It is well known that the key differences between the two theories can be traced back to 
the fact that dynamics of LQC incorporates the quantum nature of geometry through the 
area gap A £pj. The question then is: Can one regard the WDW theory as the limit of LQC 
when quantum geometry effects are ignored by taking the mathematical limit A — > 0? Since 
it is analytically soluble, sLQC is well suited to probe this issue. We used the complete set 
of Dirac observables, p^\, to compare the two theories and obtained two results. First, 
suppose we fix any semi- infinite interval / of internal time and an e > 0. Then one can 
approximate sLQC by the WDW theory to within e over the time interval interval I simply 
by shrinking by hand the area gap sufficiently. In this sense the answer to the question is in 
the affirmative. However, if one is interested in the global time evolution — i.e., if we let / be 
the full real line — then the answer is in the negative: No matter how much we shrink the area 
gap, if we wait long enough the difference in the predictions of the two theories will become 
as large as we want. Furthermore, in terms of global behavior in time, sLQC fails to admit 
any well-defined limit as the area gap shrinks to zero. In this sense sLQC is a fundamentally 
discrete theory. The use of a non-zero A in LQC is not just an intermediate step to make 
the quantum theory mathematically manageable. Rather, the presence of a non-zero area 
gap is a central physical feature of LQC, an imprint of the quantum geometry of full LQG on 
this symmetry reduced theory. The dramatic difference from the more familiar lattice gauge 
theories is that whereas QCD can be meaningfully formulated in the continuum, quantum 
general relativity cannot; in LQG it is essential to use quantum geometry at the Planck 
scale. 

We wish to emphasize, however, that the analysis of this paper has two main limitations. 
First, we restricted ourselves to the simplest model, the k=0 FRW cosmology coupled to 
a massless scalar field. The robustness refers to states within this model (although most 
results can be probably generalized by including a cosmological constant or allowing the 
k—1 closed models). The second and much more important limitation is that as of now 
LQC has not been systematically derived from LQG. In the LQC Hamiltonian constraint, 
the area gap enters through the definition of the field strength F % ah via holonomies of the 
gravitational connection A l a . This strategy is standard in full LQG. However, in LQC one 
breaks diffeomorphism invariance through gauge fixing. While this is the standard practice 
in all approaches to (classical and) quantum cosmology, it implies that we cannot directly 
employ the full strategy used so far in LQG. The new element — which may also be useful 
in a suitably gauge fixed version of full LQG — is that the loop along which the holonomy 
is defined is shrunk not to zero but only till it encloses an area a Q of Planck size. Now the 
most natural value of a Q is the smallest non-zero area eigenvalue A £p Y in the class of states 
relevant to LQC, called the area gap. However, in this step one parachutes by hand a result 
from full LQG into LQC. This is analogous to the procedure used in the Bohr atom where 
one puts in the quantization j = nh of angular momentum by hand. In retrospect, this step 
is parachuted from the more complete quantum mechanics. However, a fuller understanding 
shows that while angular momentum is indeed quantized, its eigenvalues are \/ j(j + 1) h 
rather than nh. In a similar vein the quantization of area is parachuted into LQC from 
the more complete theory of quantum geometry in LQG. We believe that this step will be 
eventually justified through a systematic derivation of LQC from LQG. However, just as the 
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correct eigenvalues of angular momentum are more subtle than the 'natural' or 'obvious' 
values nh used in the Bohr atom, the correct value of the area a a to which the holonomy 
loops have to be shrunk may have a value different from A£p V For instance, arriving at 
LQC from LQG may well require a coarse graining which could lead to a 'dressing' of A. 

We will conclude by comparing and contrasting our model with another soluble model 
that has recently appeared in the literature also for the k = cosmology with a massless 
scalar field [15]. There, in essence one begins with a classical constraint, c 2 p 2 = const pf^, 
takes its positive square root and quantizes by replacing c by sine and p by —ihd c . Thus 
one first simplifies the classical constraint and then quantizes, using rules motivated by 
LQC. 11 By contrast, in this paper started with the lapse function suited to harmonic rather 
than proper time but then followed the same procedure as in [3] to obtain the improved 
quantum constraint. Therefore, reservations [15] that the simplified model may leave out 
interesting physics are inapplicable for the model discussed in this paper. At least for the k=0 
FRW models with a massless scalar field, the bounce scenario that emerged from numerical 
simulations in [3] appears to be robust, i.e., is not tied to states which are semi-classical at 
late times. 

Finally, there are three other key differences between the two exactly soluble models. 
First, the simplified constraint used in [15] corresponds to a difference equation with uni- 
form steps in p ~ a 2 rather than in volume (~ a 3 ). Therefore, it is the analog of the older 
Hamiltonian constraint [10, 11] used in LQC which, as we discussed in sections I and III A, 
leads to physically unacceptable predictions. Second, simplifications of [15] removes, by 
hand, the right or the left moving sector of the theory. This truncation is motivated by con- 
siderations of mathematical simplicity but does not appear to have any physical justification. 
Third, as in LQC, the 'evolution equation' satisfied by physical states in our simplified model 
involves the square-root of a positive operator (0) and is thus fundamentally non-local. This 
is a direct consequence of the group averaging procedure used to obtain the physical Hilbert 
space. By contrast, because the square-root of the constraint is taken classically in [15], the 
resulting operator is local (in c) unlike the corresponding operator in full LQC [3]. In the 
terminology of section III, while the physical states of the model discussed in this paper 
have only 'positive frequency' modes but with both right and left moving components, those 
of [15] contain both positive and negative frequency modes but restricted to be either left 
or right moving. 

The simplified model introduced in this paper will be used in [19] also to analyze proper- 
ties of fluctuations in volume. Again it will turn out that indications provided by numerical 
simulations involving states which are semi-classical at late times [3] are realized by generic 
states. 



1 The replacement c — s- sine is motivated by the older, 'fi quantization' of LQC [2, 10, 11] but with /i 
set to 1. However, since it is just c rather than c 2 that now appears in the square root of the classical 
constraint, it cannot be directly related to the field strength F* b , whence the replacement c — > sine no 
longer 'descends' from LQG. See the remark at the end of section III. 
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Appendix A: Starting with harmonic versus proper time already in the classical 
theory 

In the 'improved' dynamics of [3], one began with proper time in the classical theory 
(where the lapse is given by N = 1) and switched to the scalar field time only after obtaining 
the quantum Hamiltonian constraint. In section III by contrast we first adapted the classical 
theory to the scalar field time (by using N = a 3 ) and then proceeded with quantization. The 
two strategies yield slightly different factor orderings in the quantum Hamiltonian constraint. 
In this Appendix we will discuss the detailed relation between them. It will turn out that 
the difference is negligible on physical grounds. 

If one begins with the lapse N = 1 and then passes to quantum theory, the gravitational 
part of the 'improved' Hamiltonian is given by (see Eqs. (2.23)-(2.25) of [3]): 

Asmv , <j>) = (sin Ab) A(u) (sin Ab) ^/(u, <j>) (Al) 

with 

A(u) = -^-^\u\ \\ v + \\-\ v -\\ , (A2) 

7A' 3 1 1 

where A > is given by A 2 = A £ P1 = 2-^/3^7 ^pi- The matter part of the Hamiltonian 
constraint is given by: 

CUtt <f>)) = 8nGp 2 w V(v, <t>) ■ (A3) 

Using the Thiemann strategy, [25, 26], we can calculate the operator V~ Y in the v- 
representation. As expected, this operator is diagonal: 

- — - ~ 27 1 2 2 ,-i ~ 

V-^(u) = — — - |^ + A|3 - \ u - A|3 3 ^(z/,0) 
64 A 6 a 1 1 

=: B(uMu,c/>). (A4) 

Thus, if one begins with N = 1 in the classical theory, the Hamiltonian constraint is given 
by 

^^(i/) sinAbA(z/) sinAb^M) (A5) 

The key difference from the constraint (3.7) used in the main text is that in (3.7) each of 
the functions A(v) and (i?(^))~ 1 is replaced by (certain multiples of) \u\. 
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The question then is: Under what conditions does (A5) reduce to (3.7)? Note first that 
in the main text we restricted ourselves to the lattice v = 4n\ = 4ny/A£p\. On the points 
of this lattice, we have A{y) = —(12ir£p l / , yX 2 ) |z/|. For B(v) let us make a weak assumption: 

• Let us replace B(u) with its WDW value [2-n^i\ x \v\)~ x . This amounts to assuming 
0{\\/u\) <C 1. Now, because of the form of the inner product, the state \u) = on 
which the approximation would have been the worst does not feature in the physical 
Hilbert space of [3] nor in the physical Hilbert space used in the main text. As figure 1 
shows the 'error' is only 1.43% for v = 4A, 0.02% for v = 8A and decreases extremely 
rapidly for higher v. 
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B(v)o.4 
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0.2 
0.1 




0* ' ' ' ' ' 

2 4 6 8 10 

v/X 



FIG. 1: Crosses denote values of the LQC function B{u) for u/\ = n. The continuous curve 
represents the sLQC approximation used in this paper. Physically relevant points are v/\ = An 
with n > 1. The relative 'error' is 1.43% for n = 1, 0.02% for n = 2 and further decreases extremely 
rapidly for higher n. 

With this simplification 12 , the total constraint (C^rav + C'matt) 1 ^ = takes the form 

02 t i i\ o ^ i i sinAb sinAb ~ JN 

= 37rC |^| — ^- — |^| — ^ — ^(^, 0) (A6) 

which is precisely the constraint (3.8) used in the main text. To summarize, mathematically 
the Hamiltonian constraint of [3] reduces to the one used in the main text if one replaces 
B{y) in [3] by its WDW value. The difference between the two constraint is less than 2% 
even for the state concentrated at v = 4A and decreases extremely rapidly for states with 
support at higher values of v. 



In [29] properties of the Hamiltonian constraint of [3] were established using a procedure which began 
precisely with our simplified operator. That analysis also provides a systematic leading order correction 
to the 'simplification' made here. 
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Finally, it is natural to ask if the simplifying assumption on B(v) is violated near the 
bounce point. If the violation is significant, then the conclusions drawn from the solvable 
model could be qualitatively different from those drawn from the constraint used in [3]. 
Now, since in the v representation physical states have support only on points v = 4nA 
with \n\ > 0, it follows that V m [ n is necessarily bounded below by 87T7^pjA. This value is 
in fact attained by the physical state x( u i 0) = (expiv^©^) 5|n|,i in the //-representation, 
or, 0) = (expiy/Otfi) sin(4 arctan \JYItxG x) in the x representation. In this state, the 
expectation value {p (</,)) of the constant of motion p^ is ~ h~ (in the classical units G=c=l). 
Therefore this state belongs to the extreme quantum regime. For example, in the case of 
k=l, closed universes, classical Einstein's equations imply that a universe with such a value 
of p^) can expand out to a maximum radius only of ~ 0.3£p\ before recollapsing to a big- 
crunch. Even for such an extreme quantum state, the treatment of the main text is close to 
that in [3]: as figure 1 shows, for B(v) the relative error is ~ 1.4%. In any state which can 
be thought of as representing a classical universe at late times even in a very weak sense, 
the difference between the two factor orderings would be further suppressed by an enormous 
factor. 
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